Distance growth of quantum states due to initial system— environment correlations 
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Intriguing features of the distance between two arbitrary states of an open quantum system are 
identified that are induced by initial system-environment correlations. As an example, we analyze 
a qubit dephasingly coupled to a bosonic environment. Within tailored parameter regimes, initial 
correlations are shown to substantially increase a distance between two qubit states evolving to 
long-time limit states according to exact non-Markovian dynamics. It exemplifies the breakdown of 
the distance contractivity of the reduced dynamics. 
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Introduction. -There are various states of physical sys- 
tems changing and evolving in time. Thermodynamic 
states or quantum states are well known examples. One 
can pose natural questions: How "similar" is a given 
nonequilibrium state to the equilibrium one? How " sim- 
ilar" is an entangled quantum state to the unentangled 
one? Is there a natural quantifier to measure the sep- 
aration between states? The answer is not unique and 
depends on a character of a problem which we analyze. 
In the simplest way, we could say that two states are close 
to each other if averaged values of observables in these 
states are close to each other. For quantum systems, 
there is a natural distance induced by a norm which in 
turn is induced by the scalar product of vectors in Hilbert 
space of a given quantum system. One of the central is- 
sues of quantum engineering is to prepare a quantum 
system in a state having desired features. Doing this 
one is faced with many difficulties; unavoidable decoher- 
ence, imperfect implementation of quantum operations 
and many other sources of unpredictable mess. The sta- 
tistical similarity of quantum states can be related to cer- 
tain measures of a distance between them [l| . One of the 
most popular quantifiers for the distiguishability of the 
probability distributions arising from the measurements 
conducted on quantum states (or shortly 'distinguishabil- 
ity of states') is the trace distance Q- The trace distance 
of any two states (density matrices) pi and p2 is defined 
as 



D[pi,p2] = IIpi -P2I 



(1) 



where for any operator A its norm is defined by the re- 
lation ||yl|| = (1/2) TtVA'<A. Any positive and trace- 
preserving map £ defined on the whole space of trace 
class operators on the Hilbert space is contractive, i.e. 



D[£{pi),£{p2)]<D[p,,p2]. 



(2) 



In particular, when £ = £t is a completely positive quan- 
tum dynamical semigroup such that p{t) — £tp{Q), then 

D[p^{t),p2{t)]<D[p^{s),p2{s)] for i>s, (3) 

i.e. the distance cannot increase in time. It means that 
the distinguishability of any states can not increase above 
an initial value. Motivated by findings recently reported 



in Ref. [3|, we study an open quantum system which is 
initially correlated with its environment [J| . We consider 
a dephasing model [5|, l6| of decoherence of a qubit inter- 
acting with a bosonic environment. If the qubit and its 
environment is initially prepared in an uncorrelated state, 
the reduced dynamics is completely positive and hence 
contractive. In the presence of initial correlations one is 
left with positive (but non completely positive) dynam- 
ics [3, Q and the contractivity may fail f3| . The model 
of dephasing allows to find an exact reduced dynamics 
and reveals highly non-trivial properties of distinguisha- 
bility with respect to the trace distance. We explicitly 
show that under tailored regimes the contractivity fails 
for some initially entangled states. The trace distance of 
different states grows above its initial value even in the 
long time limit. In consequence, the distinguishability of 
the long-time limit states can increase above its initial 
value. 

Model. -The system we study is a qubit Q, formed by 
an arbitrary two-level system coupled to its environment 
B. We consider the case when the process of energy 
dissipation is negligible and only pure decoherence is a 
mechanism for decoherence of the qubit. It leads to an 
irreversible process of information loss. We model such 
a system by the Hamiltonian [h—l) 

H^Hq®Ib+Iq®Hb + S''®Hi, (4) 

HQ=eS\ Hb= dujh{uj)a\uj)a{uj), (5) 

"'0 



Hr 



doj [g*{uj)a{uj) + g{uj)a\uj)\ 



(6) 



where S^ is the z-componcnt of the spin operator and is 
represented by the diagonal matrix S^ = diag[l, —1] of 
elements 1 and —1. The parameter e is the qubit energy 
splitting, Iq and I_b are identity operators (matrices) in 
corresponding Hilbert spaces of the qubit Q and the en- 
vironment B, respectively. The operators a^(aj) and a{Lu) 
are the bosonic creation and annihilation operators, re- 
spectively. The real- valued spectrum function h{uj) char- 
acterizes the environment. The coupling is described by 
the function g(uj) and the function g*{uj) is the complex 
conjugate to g{uj). The Hamiltonian ^ can be rewritten 



in the block-diagonal structure Q, 

H = diag[H+,H-], H± = Hb ± Hi ± els- (7) 

We consider a correlated initial state of the qubit- 
environment system in the form similar to that in Ref. 
[lO|, namely, 



\^iO))=b+\l)®\no)+b-\-l)®\Qx). 



(8) 



The states |1) and | — 1) denote the excited and ground 
states of the qubit, respectively, the non-zero complex 
numbers b+ and b- are such that |6+p + |6_p = 1. The 
states l^lo) and \i^\) are environment states. As an ex- 
ample, let us assume that jfio) is an environment ground 
state and 



Na>=C,-i [(i_A)|f7o)+A|%)] 



(9) 



is a linear combination of a ground state \flo) and the 
coherent state jJl/) = D{f)\^o), where the displacement 
operator D{f) reads [11| 

D{f) = exp |y" doj [f{io)a^{oj) - r{Lo)a{u)] \ (10) 

for an arbitrary square-integrable function /. The con- 
stant C\ normalizes the state \fl\) and reads 

Cl = {l- A)2 + A^ + 2A(1 - X)Re{no\nf), (11) 

where Re is a real part of the scalar product {QolQf) 
of two states in the environment Hilbert space. The pa- 
rameter A € [0, 1] controls the initial entanglement of the 
qubit and environment. For A = the qubit and envi- 
ronment are initially uncorrelated while for A = 1 the 
entanglement is, for a given class of initial states, maxi- 
mal. The states with A = 1 are not maximally entangled 
in the usual sense as the coherent states, forming an over- 
complete set of states, are not mutually orthogonal. The 
state of the total system at time i > has the form 

I^I'(i)) = b+\l) ® |^+(i)) +6_| - 1) ® |^_(t)), (12) 

where |V'+(0) — exp{—iH+t)\ilo) and \tp^{t)) = 
e'Kjp(—iH-t)\n\). The density matrix px{t) of the qubit 
is obtained from the relation 



p,{t) = TrB{\^{t)){^{t)\} 



(13) 



and Ttb denotes the partial tracing over the environ- 
ment. Its explicit form reads 



Pxit) 



|&+|2 b+b*_Ax{t) 

bXb^Alit) |6_|2 



(14) 



Ax{t) = C^^ e-2^"*e-''« 1 - A + Xe-"""^^'^ e'^'^ , (15) 
where H 



r{t) = 4 / dugiiio) [I - cos(c^i)] , 
Jo 

D 1 r°° 

s(i) - 2 / du:gh{iu)f{u:) [1 - cos(wi)] - 7^ / dujf{iu), 



where gh{'^) — 9{'^)/h{Lu) and 



$(i) = / dujghiu})f{uj)sin{ujt). 
Jo 



(16) 



For the sake of simplicity, we have assumed that the func- 
tions g{uj) and /(w) are real. The modeling of the qubit 
coupling to the environment is performed in terms of the 
spectral density function 



.g;j(w) = aw^ exp(-a;/wc). 



(17) 



where a > is the qubit- environment coupling constant, 
/i > — 1 is the "ohmicity" parameter (the case p — 
corresponds to the ohmic and /x > to super-ohmic en- 
vironments, respectively) and Wc is a cut-off frequency. 
The coherent state \flf) is determined by the function 



f^{uj) ^ -/uj" ^exp(-w/iLJc)- 



(18) 



This choice is arbitrary but convenient as it allows for 
finding an explicte analytic formulas for 

r(t)=4£(a,^,i)], 
s{t) - 2/:(VS7, {p + i^)/2, t) - 7^(^^)c^,^ 
cos [/i arctan(wc^)] 



$(i) = ^/a^T{K)uj^ 



(1-1-^2^2)^/2 

^ sin [k arctan(a;ct)] 

(1+^2^2)^/2 



(19) 



where r(i/) is the Euler gamma function and the param- 
eter K = {p + v)/2. 

Distance of states. -Let us examine the distinguisha- 
bility between two various, initially correlated (A ^ 0) 
states. The corresponding trace distance reads 



(20) 



7?^K(i),PA.(i)] = (|6V'^|'-|&f 
+ \b^:h^l>A,At)-b^:'b^^>A,.At)? 



As the initial state ^ of the qubit-environment compos- 
ite system is generically entangled, the non-local opera- 
tions are involved in its preparation. Such a procedure 
would essentially require highly sophisticated quantum 
engineering. Naively speaking one can recognize the fol- 
lowing steps: (i) the first, when one prepares the state 
\fl\) given by Eq. ©, next "tensorize" it with the ground 
state of the qubit | — 1); (H) the second step, when one 
superpose, with the weights b±, the result of the first 
step with the excited qubit |1) state tensorized with the 
vacuum \Qo)- Both steps demand using precise technics 
and their details are essentially beyond the scope of this 
paper. A different preparations of initial states ([U can 
be provided both via changing details of the first step 
and the second step of the procedure. In other words, 
one can change the parameters b± in ([8]) or manipulate 
on the state Iflx) (by changing A). When two different 
states are determined by two different sets of numbers 
bj. {k = 1,2) in the superposition ([5]) with the same 



state \Q\] 
form 



then Axi (i) = Ax^ (i) and Eq. ^0^ takes the 



D'[pxAt),Px.{t)] = [\b^^^\' 



-16 



(2) ,2 



l^w^a)* 



k(^)b(^>l2 



\AxAt)? 



(21) 



The function |A^j(i)| is a decaying function of time and 
the distance (PT|) also decays as time grows. In conse- 
quence the distinguishabihty between two initially cor- 
related states always becomes lowered. It is clear that 

(k) 

no matter what bj_ are, the reduced dynamics results in 
contraction of the distance between the states. So, the 
only chance to observe the growth of the distinguishabil- 
ity of qubit reduced states is to modify the parameters of 
environment encoded in \flx) in Eq.®. Let us consider 
this case assuming that &j_ = b± are fixed and the same 
for k = 1,2, and two different states are determined by 
two different sets of numbers Ai and A2 . In such a case 
the distance reads 

\b+b*_\-' D[pxAt),PxAt)] = \AxAt) - AxAt)\ 
= e-'-(*) {a^ _^ ^2^2s(t) ^ 2a6e"(*) cos[2$(i)]} , (22) 



where 



1 - Ai 1 - A2 



Cx 



Cx 



b = 



Cx, 



A. 



(23) 



In Fig. [lla), we show time evolution of the trace distance 
D[px{t), pa{t)] between the initially correlated and non- 
correlated states for four selected degrees of correlation 
determined by the value of the correlation parameter A. 
At time t = 0, the distance D = 0.004,0.03,0.08,0.19 
for A = 0.1,0.25,0.4,0.6, respectively. In the long 
time limit, the distance between states takes the val- 
ues D = 0.03,0.08,0.12,0.15, respectively. In the first 
three cases, the distance between long-time limit states 
is greater than the initial distance and the distinguisha- 
bihty of final states is better than initial states. There 
is some optimal correlation A for which the distinguisha- 
bihty of final states is the best. However, if the initial 
correlation is strong enough, the distance for long time 
is smaller than at the initial time. It is important ob- 
servation that there exists a critical value Ac (depending 
on other parameters) such that for A > Ac the station- 
ary states are less distinguishable than initial states. In 
some regimes (as e.g. the case A = 0.25, 0.4 in Fig. [I^a)), 
at early stage of evolution, the distance decreases reach- 
ing a minimum and next it increases and saturates for 
long time. In Fig. [IJb) and Fig. [Ijc) we show how 
the distance can depend on environment characteristics; 
the coupling constant a and the " super-ohmicity" /i > 0. 
It follows that strong coupling of the qubit to environ- 
ment diminishes the distinguishabihty of initial states. 
Similarly, if the environment is more super-ohmic, the 
distinguishabihty of states is weaker. In Fig. [21 the in- 
fluence of parameters characterizing initial environment 
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FIG. 1: Time evolution of the trace distance D[px{t), po{t)] 
between the initially correlated and non-correlated states for 
selected values of parameters. Dimensionless time is in unit 
of (jjc. The remaining parameters are: e = 1, {a)p = 0.01, 
u = 0.05, a = 0.01,7 = 0.05; (b) p = 0.01, u = 0.05 , 
7 = 0.05. A = 0.25; (c) v = 0.05 , 7 = 0.05. A = 0.25, 
a — 0.01. The horizontal lines in (b) and (c) indicate the 
initial distance. 



state is depicted. Both regions, where the "distinguisha- 
bihty gain" (i. e. the distance in the long time limit is 
greater than for the initial time) is reached, are limited 
and located in the corner of the corresponding parameter 
planes. It means that not all but specific environments 
can induce the " distinguishabihty gain" . 

Finally, we consider the case when both Ai and A2 
are non-zero, see Fig. 3. At time i = 0, the distance is 
D = 0.015, 0.014, 0.0071, 0.026 for A = 0, 0.05, 0.3, 0.4, re- 
spectively. In the long time limit, the distance takes the 
values D = 0.039,0.031,0.0072,0.020, respectively. In 




FIG. 2: In the region limited by the curves and the 
horizontal axes (denoted by +), limt^oo D[px-^{t),po{t)] > 
D[pai (0),po(0)]. Outside this region (denoted by — ), the 
distance between two states at time f = is greater than 
between the long-time limit states (the contractive region). 
The parameters are: e = 1, 7 = 0.05, p — 0.01; Inset: 
/i = 0.01, i/ = 0.05. 
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FIG. 3: Time evolution of the trace distance D[pxj^{t), px^it)] 
for fixed Ai = 0.25 and several values of A2. Time is scaled 
as in Fig. 1. The parameters p — 0.01, u — 0.05 a = 0.01, 
7 = 0.25 and b+ = b- = 1/^2. 



the first three cases, the final distance is greater than the 
initial distance and the distinguishability of final states 
is better than initial states. The region of the (Ai, A2)- 
plane, where this effect occurs is bounded by Ai and A2 
axes and a decreasing function A2 vs Ai. From the pre- 
sented example we detect that if A2 exceeds some fixed 
value of Ai, this effect is distinctly absent and the dis- 
tance in the long time limit can never be greater than for 

Summary. -Recent investigations have often shown 
that dynamics of quantum systems can exhibit various 
counter-intuitive and non-trivial features which can be 
inferred from detailed analysis of highly sophisticated 
theoretical models. It is of great importance to verify 
the predictions in a carefully prepared experiments. An 
effective design of such experiments would require guide- 
lines provided by theoretical studies of realistic physical 
models. Our work reported in this paper, as a step in 
this direction, could be placed somewhere in between: 
we consider a fairly simple model of decoherence under 
very specific conditions but present results which are ex- 
act. Despite its simplicity the model is realistic enough to 
be experimentally accessible \^. According to common 
wisdom the decoherence resulting from the "openness" of 
quantum systems causes blurring of the information en- 
coded in quantum states: as it has recently been shown 
[3| in the context of distinguishability of quantum states, 
it is not always the case. In this paper we provide an 
explicit realization of the recent suggestion on the dis- 
tinguishability growth due to initial qubit-environment 
entanglement. The distinguishability growth occurs not 
only at the short time scales but is shown to be a fea- 
ture of long-time limit states. This feature seems to be 
favorable for the potential experimental verification of 
the predictions presented in \3\ especially when the de- 
sired short-time growth would occur in ultra-short time 
scales. We have shown that the induced distinguishabil- 
ity growth is not generic for the considered model and we 
identified the parameter regimes where the effect is the 
most apparent. 
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